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Введение. Разработаны математические модели для про- 
гнозирования показателей качества кокса М5 и Мио. Рас- 
четы выполнены индивидуально для каждой из коксовых 
батарей коксохимического производства (КХП) Магнито 
горского металлургического комбината (ММК). В основе 
моделирования — параметры шихты: сумма отощающих 
компонентов ОК, %; отражательная способность витрини- 
та Ко, %. Такие модели необходимы при управлении каче- 
ством продукции и оптимизации, направленной на сниже- 
ние производственных затрат. Цель исследования — по- 
строение адекватных математических моделей для про- 
гнозирования показателей качества кокса М5 и Ми в 
условиях КХП ММК. Таким образом, предполагается, что 
ММК получит собственные модели, не уступающие по 
точности прогнозирования аналогам, используемым дру- 
гими коксохимическими предприятиями России. 
Материалы и методы. В качестве универсального ап- 
проксиматора для построения математических моделей 
использовались нейронные сети. При выборе их архитек- 
туры исходили из минимального количества нейронов и 
слоев сети. Кроме того, учитывалась минимизация ошиб- 
ки прогнозирования на новой выборке, которая не исполь- 
зовалась при обучении и тестировании. 

Результаты исследования. Разработка основана на петро- 
графических параметрах шихты: суммы отощающих ком- 
понентов ОК (по ГОСТ 12112); показатель отражения 
витринита в шихте Ку (ГОСТ 12113). С помощью искус- 
ственных нейронных сетей построены одномерные мате- 
матические модели для прогнозирования показателей 
качества кокса на дробимость М>5 и истираемость Муо 
(ГОСТ 5953). Созданные модели представлены в графиче- 
ском виде. Оценены их прогнозирующие способности. 
Обсуждение и заключения. В моделях, разработанных в 
рамках данного исследования, использованы только пет- 
рографические параметры шихты. Не учитывались сово- 
купные данные по техническому и пластометрическому 
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анализам. Этим подход, представленный в данной статье, 
отличается от моделей, реализованных для других КХП, 
— например, на Нижнетагильском металлургическом 
комбинате (НТМК), Новокузнецком металлургическом 
комбинате (НКМК), Западно-Сибирском металлургиче- 
ском комбинат (ЗСМК). Несмотря на это, доказана адек- 
ватность полученных зависимостей. Предполагается ис- 
пользовать их для оптимизации петрографических пара- 
метров шихты по различным критериям оптимальности 
показателей качества кокса. 
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Введение. В настоящее время на коксохимических производствах (КХП) проводятся мероприятия по 
снижению себестоимости производимой продукции [1-3]. Оптимальные параметры этой работы на КХП ПАО 
«Магнитогорский металлургический комбинат» (ММК) могут быть определены с помощью адекватной мате- 
матической модели. Предполагается, что она будет описывать процессы превращения угольной шихты в кокс и 
формирования его прочностных свойств. 

Материалы и методы. Выбор и обоснование независимых параметров для построения математи- 
ческих моделей. В [4] на основе построения и применения адаптивной модели установлен следующий факт. На 
точность прогнозирования показателей качества кокса влияют только петрографические параметры шихты, а не 
данные технического и пластометрического анализов. В связи с этим при построении моделей для прогнозиро- 
вания показателей качества кокса М5 и Ми в качестве независимых переменных выбраны следующие петро- 
графические характеристики шихты: сумма отощающих компонентов (ОК %) и показатель отражения витрини- 
та (Ко %). 

Многолетние данные о работе КХП ММК для 7-бис, 8-бис, 13-й и14-й батарей позволяют утверждать, 
что существует тесная взаимосвязь между показателями качества кокса М›5, Мии ОК, К\ соответственно. 

Выбор и обоснование вида математического уравнения для построения модели. В работах [З, 5, 6, 7] 
было установлено, что линейные уравнения в качестве моделей неадекватно описывают зависимости прочностных 
свойств кокса от шихтовых параметров. Поэтому математическую модель будем строить в виде: 

М»5 = Г(х), (1) 
Мто = /(%), (2) 
где / — нелинейная функция; х — петрографический параметр ОК или Кб. 

В работах [3, 6-13] с помощью множественного регрессионного анализа были найдены зависимости по типу 
(1) и (2) от нескольких независимых переменных. В настоящее время искусственная нейронная сеть является самым 
современным и удобным алгоритмом обработки информации (с учетом имеющихся прикладных пакетов программ для 
ЭВМ) [14, 15]. Это объясняет применение нейронных сетей в данной работе. В промышленных условиях работы кок- 
совой батареи одному входному параметру, например ОК, может соответствовать несколько значений М5, Ми. Это 
объясняется возможной неравномерностью состава шихты, насыпной плотности по ширине и высоте камеры коксова- 
ния и другими факторами. В таком случае нейронные сети будут выступать в качестве сглаживателя данных [15]. Под 
термином «искусственная нейронная сеть» («нейронная сеть») в представленной работе будет подразумеваться зависи- 
мость по типу (1) или (2). 

Результаты исследования. Построение математических моделей, проверка их адекватности и работо- 
способности. Данные КХП ММК фиксировались в течение 2011-2013 годов. На их основе с помощью программы 
МАТГАВ в пакете М№еига! Мегрогк Тоо! для каждой коксовой батареи (табл. 1) подбиралась своя полносвязная нейронная 
сеть. В такой сети отсутствуют пропуски связи между слоями нейронов и сигнал (числовое значение) с каждого нейро- 
на из предыдущего слоя участвует в формировании сигнала (числового значения) каждого нейрона в следующем слое. 
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Таким образом, следует подчеркнуть, что при создании моделей для коксовых батарей не использовался инструмента- 
рий встраивания априорной информации об объекте исследования в нейронные сети. Они применялись для современ- 
ного и удобного нахождения зависимостей между независимыми параметрами и показателями качества кокса. 
Таблица 1 
Та 1 
Характеристики подобранных одномерных математических моделей для прогнозирования показателей 


качества кокса Мь5, Мо для коксовых батарей, обслуживаемых 1-м блоком углеподготовительного центра 
(УПЦ) КХП ММК 


Срагасчет5йсу орстешей 12 тайетайса! то4е юг рге@сипе М»; апа Мтосоке диа Шу таГсеу [ог соке Байетеху 
орегшеа Бу ше 151 Боск ор соа[ ргераганоп сетте (СРС) оГССР ММК 


Номер коксовой батареи 


Характеристика 
М 


Независимая переменная, от 
которой строилась модель 


д Е ь 
сети 
С 6 
реднее а ошибки 0.0275 0,0355 | —0,0046 | -0,0215 0,0146 0,0093 
описания 
С - 
реднеквадратическое откло 0.5954 0.3952 0,6212 0,3910 0,9311 0,3159 
нение ошибки описания 


Вероятность того, что мера 0,00001 (14) 
расхождения /° для ошибки 0.2312 0,00000 (10) 
описания не меньше, чем | 0,3063 (12) 0,00001 (20) 0,5175 (14) 
зафиксированная” о 0,00003 (50) 
0,00001 (100)* 


Интервал (0.61: 0,63) 
попадания 
ошибки опи- 
а (-0,94; 0,96) 
Значение ошибки прогнози- 
рования на новой выборке 
(месячные данные за 2014 
год): минимальное отрица- 
тельное / максимальное по- 


ложительное 





`В ошибки описания с нормальным законом распределения [16].обозначении архитектуры нейронной сети литера { обозначает, что в каче- 
стве передаточной функции на данном слое использовался гиперболический тангенс. Например, архитектура 11—3—И означает, что в сети 
имеется три слоя нейронов: {1-— входной слой, {3 скрытый слой, — выходной слой. Во входном слое имеется 1 нейрон, в скрытом слое — 3 
нейрона, в выходном — 1. Для прогнозирования параметра М>5 для батареи 8-бис была подобрана особенная архитектура: в скрытом слое 
не три, как у всех остальных моделей, а два нейрона. 

* Ошибка описания — принятый в данной статье термин, обозначающий разницу между прогнозируемым значением какого-либо показате- 
ля качества кокса, полученным с помощью математической модели, и фактическим из выборки, на которой строилась модель. 

? В скобках указано число разрядов, на котором строился статистический ряд. 

* При данном значении вероятности можно утверждать, что распределение ошибки описания не согласуется с нормальным законом распре- 
деления. Значение приводимой вероятности может несколько зависеть от количества разрядов, на котором строился статистический ряд. 
Поэтому приведено несколько значений вероятности в зависимости от числа разрядов, для того чтобы в полной мере убедиться в адекват- 
ности сделанного вывода о несогласованности тестового множества, которое было случайным образом сформировано из входных данных [15]. Для 
каждой модели была выбрана только одна независимая переменная (ОК или Ко), и нейронные сети использовались для сглаживания данных. Учитывая 
это, нейронная сеть заведомо не должна быть большой. Для малых нейронных сетей наиболее эффективен (сходимость, затраты машинного времени) 
алгоритм обучения Левенберга — Марквардта [15], который и был использован в данной работе при построении моделей. Поскольку результаты обу- 
чения зависят от начальных значений весов, то обучение на одной архитектуре нейронной сети происходило несколько раз после случайной инициали- 
зации весов. В качестве конечного результата оставлялись только те наборы коэффициентов в нейронной сети, при которых для данной архитектуры 
было наименьшее значение дисперсии (теап 5диаге етот — МЕ в МАТЕАВ) в обучающей выборке. Коэффициент обучения не изменялся, использова- 
лись настройки по умолчанию. Всего с целью обучения и нахождения нейронных сетей было взято по 157—168 значений для каждой батареи. 
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Для выбора нейронной сети использовались следующие дополнительные критерии: 

— подтверждение гипотезы о соответствии величины ошибки описания нормальному закону распределения по крите- 
риюх; 

— попадание максимального по модулю значения ошибки прогнозирования” на новой выборке (месячные данные 
КХП ММК за 2014 год) в интервал ошибки описания для данной модели (см. табл. 1). 

В ходе поиска модели для 13—14-Й батареи по показателю М»5 не удалось подобрать соответствующую 
нейронную сеть, чтобы ошибка описания соответствовала нормальному распределению. Однако применение данной 
модели к новой выборке (данные 2014 года) показало удовлетворительную прогнозирующую способность. Отмечен- 
ные значения аппроксимированы с помощью нейронных сетей. Это позволило прогнозировать показатели качества 
кокса М›5 и Мио для коксовых батарей 7-бис, 8-бис, 13-Й и 14-й батарей. Результаты приведены на рис. 1-4. 
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Рис. 1. Результаты аппроксимации отмеченных значений показателя дробимости кокса М>5 по сумме отощающих 
компонентов шихты ОК для 7-бис и 8-бис батарей с помощью нейронных сетей 
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8,60 
8,50 
8,40 
М» % 8,30 
8,20 
8,10 


8,00 





ОК, % 
Рис. 2. Результаты аппроксимации отмеченных значений показателя истираемости кокса Мио по сумме отощающих 
компонентов шихты ОК для 7-бис и 8-бис батарей с помощью нейронных сетей 


Ето. 2. Арргомтайоп гезий5 оГоб5егуей уамез ор Муту соке абтазтоп таех Бу 5ит оГтет сотропеп оГОК сйагее 
юг 7-65 апа 89-615 Байетеу изтз пеига[ пелот/$ 
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Рис. 3. Результаты аппроксимации отмеченных значений показателя дробимости М» кокса по показателю отражения 
витринита шихты А‹ для 13—14-й батарей с помощью нейронных сетей 
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"Под термином «ошибка прогнозирования» в данной работе подразумевается разница между прогнозируемым значением какого-либо пока- 
зателя качества кокса, полученным с помощью математической модели, и фактическим из новой выборки, которая не использовалась для 
построения модели. 
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Рис. 4. Результаты аппроксимации наблюденных значений для показателя истираемости Мио кокса от показателя отраже- 


ния витринита шихты К. для 13—14-й батареи с помощью нейронных сетей 
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Итак, для 7-бис и 8-бис батарей (рис. 1, 2) видно, что с ростом ОК значение М»; уменьшается, а Ми 
увеличивается. Это соответствует методике определения М›5 и Мио по ГОСТ 5953-93 (ИСО 556-80). Аналогич- 
ные рассуждения можно провести для 13-—14-й батарей (рис. 3, 4). Таким образом, модели, представленные на 
рис. 1-4, можно считать адекватными, исходя из их непротиворечия описанию зависимых переменных, прогно- 
зированию на новой выборке и в сравнении (форма кривых) с моделями других коксохимических производств 


[11, 12, 17]. 


Обсуждение и заключения. Сравнение прогнозирующих способностей полученных математиче- 
ских моделей с аналогичными моделями, построенными для других коксохимических производств. В 
табл. 2 приведены прогнозирующие способности некоторых математических моделей, относящихся к разным 
КХП, исходя из их ошибок описания и прогнозирования. 


Таблица 2 
ТаБЕ 2 


Прогнозирующая способность некоторых математических моделей для показателей качества кокса, 
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С учетом тестирования на новой выборке по данным за 2014 год прогнозирующая способность моде- 


° В оригинальной работе подобная характеристика не приводится. При расчете исходили из удвоенного значения среднеквадратического 


отклонения (правило 


двух сигм). 


7В скобках указаны остаточные дисперсии, которые приводились в оригинальной работе. 
$ ООО «ЕвразХолдинг» — компания, владеющая НТМК, НКМК, ЗСМК. 
ы Приводятся характеристики для модели с лучшей прогнозирующей способностью. 
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лей, созданных для КХИ ММК по показателям М>; и Мо (табл. 1), не уступает аналогам (табл. 2). Таким обра- 
зом, модели, представленные на рис. 1—4, адекватны для нахождения оптимальных петрографических парамет- 
ров шихты 1-го блока УПЦ и прогнозирования. 

Выводы. 

1. Показана возможность построения математических моделей на основе только петрографических 
параметров шихты. 

2. Результаты прогнозирования оценивались на выборке, которая не использовалась для построения 
моделей, и сравнивались с аналогичными моделями для других КХП [7-13]. Таким образом установлена адек- 
ватность построенных математических моделей (рис. 1-4). 

3. Построенные математические модели для прогнозирования показателей качества кокса М5; и 
Мио (рис. 1-4) адекватны для проведения оптимизации [6—10, 12, 13]. 
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